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Pbeface. 

Part II. of the present paper was written, very much in the form in which it is 
now presented, in the summer of 1913, and began with the remark in § 11, which 
appears to disprove a statement made by Poincare in regard to the convergence of 
an astronomical series. It was laid aside partly because a good deal of the work is 
only of the nature of elementary algebra, partly because the matrix notation 
employed does not seem to find favour in its application to differential equations. 
Various circumstances have, however, led me to take up the matter again, and my 
original conviction that the method of Part II. is of importance has been strengthened 
by comparing it with the less formal methods which, for the sake of introducing the 
subject, I have followed in Part I. I hope, therefore, that the following exposition 
may be thought worth while. Part III. has only the value of a concluding remark. 

The table of contents above may serve to give an idea of the scope and arrangement 
of the paper. 

PAET I. 

§ 1. Consider a linear differential equation 

U^+V^+WX-0, 

dr dr 

where U, V, W are power series in a small quantity, X, of the forms 

V = v-V\v^ +X\^2 +-^.> 
in which each of u^, v^, w^ is a linear function of 



'Vi ^ri 



P-r ?r~2 Pr-4: P4-r P2-r P- 



■r 

3 



^ denoting e\ Thus each of u.^n, v^n, W2n will contain a term independent of ^ ; we 
speak of these as the absolute terms. It is important that W contains no term 
in X^ ; and it is assumed that the quantity vju, which is independent of f , is not' a 
positive or negative integer, and that ii, v are not both zero. 

We prove that if the absolute terms in W, that is the absolute terms in 

102', w^^ tv^^ . . , , 

be suitably determined, the differential equation possesses a solution of the form 

Jv = 1 + \(pi + X (p2~^ '" ^ 
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wherein cj>^ is a linear function of f"*, f*'"^ f*'^^, ..., f^"**, f^""*, f "*", and this is a periodic 
solution. Its period is 2Tri ; we can, however, if we wish, express the same result 
with a period 2it by writing r = it. 

For the substitution of the assumed form for X requires the identity 

which, equating the coefficient of X"" to zero, will be true if 

In particular for n ~ 1 

If herein we suppose 

01 = Aif+A_if-\ tc;, = Cl^f c_if-\ 

t^, t?, Ci, c_i being given constants, we obtain 

u (Ai^+ A_iro +^ (Aif- A.^r^) +ci^+c_ir^ = o, 



which is satisfied by 



xjL| — ""* ————— J ^^^ ^ ^ 



c 



1 



For n = 2 the condition is 
Writing 

and assuming a form 
the condition becomes 

equating the coefficients of ^^, ^~^, f^ to zero, we obtain 

{4:U-\-2v) ^2= — aiAi — 6iAi— CiAi— C2, 
(4'u-2t;) A„2 = -a_iA_i + 6_iA_i— c_iA_i-c_2, 

C2 = — aiA^i-~a_iAi + 6iA_i— &_iAi— CiA_i— c_iAi, 

T 2 
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which, as vjii is not 2 or —2, determine A^, A^^^^? and G2, the last being expressible by 
means of the given coefficients of u, 1;, Ui^ Vi, tVj. 

Proceeding similarly with the general value of n, y/e at once reach the conclusion 
stated, the absolute term in tv^^ being determined in terms of the coefficients in 

§ 2. Now consider an equation 

(It'^ d'T 

where, with ^ =^ 6"", A, B, C have the forms 

which are periodic functions of t, with period 27ri^ capable of being arranged as 
power series in a parameter X, the coefficient of X'' being a linear function of 

^r Pr-2 P2-r P-r 

In accordance with the well-known theory of linear differential equations with 

a/ -i 

periodic coefficients, we substitute 



X ^ e"'" 



where /c is a constant, and so obtain a diffi3rential equation 

AX'' + 2 (/cA + B) X' + (A// + 2B/C + C) X = 0, 

which, when k is properly chosen, is to be satisfied by a periodic function X. That 
this is possible follows at once from § 1, as we now explain. 

First we can draw some inference as to the form of /c. For compare the original 
differential equation in x with the equation obtained from it by changing the sign of 
X in each of the series A, B, C It is clear that the new differential equation may 
equally be obtained from the original equation by change of r into t + tt^*, which 
changes f into ~-^; this latter change, however, only multiplies the factor e''^ by the 
constant e^'''' ; the factors e"^ appropriate to the two independent solutions of the new 
.differential equation are thus the sa-me, in their aggregate, as the factors for the 
original equation. Thus the change of the sign of X changes the two factors e"'' 
appropriate to the two independent solutions of the original differential equation 
among themselves, either by leaving both unaltered or by interchanging them. 
Assuming that /c is capable of expression as a power series in X, 

K = Kq-J- /CjX + ^gX + . . . , 
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the case in which each k is unaltered by change of the sign of X is the case in which 
only even powers enter in this series. The case in which the two values of k are 
interchanged by change of the sign of X may arise when the differential equation is 
such that for X = the two values of k are equal or differ by an integer ; in this case 
^«'Y^«T |g ^ periodic function for X = 0, and the factors e"'"", e""" do not individualise the 
functions with which they are associated. 

In the present case, the equation reduces when X = 0, to 

<^o -T-^ + 26o -7- + CqX = 0, 

dr dr 

which, if ^0 is not zero, has the two solutions e*^^, 6""^ where a-, a have the values 
Thus if we suppose not only that % is other than zero, but also that 

is not zero or a positive or negative integer, we can assume 

Then putting 

X = 1 + X^i ~f X^^2 "!"•••? 

where ^^ is a linear function of f**, f**""^, .... i^~% f""*, the differential equation for X 
can be compared with that of § 1 . In the present case there is an unknown 
quantity k entering into the coefficient A/c + B of dX/dr, but it will be seen that in 
the equations obtained by taking the successive powers of X, each unknown coefficient 
in K in this A/c + B is determined at an earlier stage as entering in the coefficient 
A/c^ + 2B/c+0, and so enters as a known coefficient. We have 

A/c + B = [% + X(aif+a„if"^) + X^(a2fHa„2f"^+a2o) + ...] [^+'^2^^+"-] 

+ X^ [or (a^f + a^i~-' + a^o) + a^K^^ b^' + 6„2^ ~' + Kl 

■ i ... J 
and similarly, 

+ \'[<t' {a,^^ + a_,C-^+a2o) + 2<T {b,^'+b_,^-' + b,o) + 2K,{a,a + b,) 

"1 • • • , 



40* 
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the absolute term in the coefficient of X* in this being 

Thus, as in § I5 we first put 
assuming, as in 1 1 it was assumed that v/u is not an integer^ that 

is not zero or integral; then the absolute term in the coefficient of X^ determines 
/c2(aoa-+?^o)? 9.nd hence /cg, and the absolute term in the coefficient of X* similarly 
determines /C4. 

The excepted case in which /c contains odd as well as even powers of X we may 
leave aside at present. 

1 3. We may apply the preceding to the much discussed^ equation 

4.(crH2X^i cos + 2X^^3 COS 20+ *..)a; = 0. 

When X = we have the two factors c^'*^^, e~''^''\ and the general case is that in 
which e^^'^^ has not the period, 27r, of the coefficients in the differential equation, 
that is, when 2cr is not an integer. First assume this to be so. Then writing 

we obtain 

X'' + 2fVX'+ {a^-K^+2Xh^ cos + 2X^/^2 cos 20+ ...) X = 0. 

Herein assume 

/c = cr + K2k^ + /C4X* + . . . , X = 1 + X^x + X 02 + « « • J 

* For this differential equation the following list of references may be useful, though it is probably far 
from complete :-~-Mathieu, ^ Louville's J.,' XIIL (1868), p, 137; Hn.L, * Coll. Math. Works/ I. , p. 255 
(' Acta Math.,' VIII. (1886)); Adams, 'Coll. Scientific Papers/ I, p. 186, II., pp. 65, 86; Tisserand, 
* Mec. Cel.,' t. III., Ch. I. ; PoiNCAiiE, ' Meth. Nouv.,' t. IL, Ch. XYIL ; Forsyth, * Linear Differential 
Equations' (1902), p. 431; Rayleigh, * Scientific Papers,' vol. Ill, (1902), p. 1; Lindemann, ^Math. 
Annal,' Bd. XXII. (1883), p. 117; Linbstedt, ^ Astr. Nachr.,' 2503 (1883); Linbstedt, *M4moires 
de FAoad. de St.. Petersbourg,' t. XXI, No. 4; Brums, ^ Astr. Nachr,,' 2533, 2553 (1883); 
Oallandreau, 'Astr. Nachr.,' 2547 (1883); Oallandreau, 'Ann. Observ.,' Paris, XXII. (1896); 
TissERANB, 'Bull. Astr.,' t. IX. (1892); Stieltjes, 'Astr. Nachr.,' 2602, 2609 (1884); Harzer, 
'Astr. Nachr.,' 2850, 2851 (1888); Moitlton and Macmillan, 'Amer. J.,' XXXIII. (1911-); 
MOULTON, 'Rendic. Palermo,' XXXII. (1911); Moulton, 'Math. Ann.,' LXXIII. (1913); Whittaker, 
'.Cambridge Congress' (1912), I, p. 366; Whittaker, Young and Milne, 'Edinburgh Math. Soc.,' 
XXXIL, 1913™14; Inge, 'Monthly Not.,' Roy. Astr. Soc, LXXV. (1915); Poincare, 'Bull, Astr.,' 
XVII. (1900). 
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where cf>^ is an integral polynomial of order r in ^ and ^~\ the quantity e^^ being 
denoted by f. Then we have 



n 



i:X''<p"„ + 2i{o-+ K^X^ + . . . ) zx>', 

+ [-2(r/c2X'-(2o-/c4 + ,c/) X''+... + 2X^1 COS 0+...][l + 2A>J = 0. 
The terms in X give 

^^+2^v.^'i+^i(f+ri) = o, 

Avhich, if we denote {cr + rf—a-^ or r{2a-+r) by u^, so that the result of substituting 
f for in <p" + 2iG-(p' is —Urf, leads to 

The terms in X^ give 

which, if we write 
leads to 

and 



2(7 \«i M_i./ o-(4o-^— l) 

By the terms in X^, X"*, we similarly find the coefficients in 

,/,3 = A3^H A_,r«+B,^+B_,rs 



and also 

60a-^--35a-' + 2 , , . 3 



• 2 



J 



Ka — 



4a-Ha-^-l)(4c7^-l)'^^''"^2a-(a-^-l)(4a-^~-l)*^'^^ 4^(^2^l)^/- 



If we change the notation, writing Q — 2t, 2ar = n, so that the differential equation 
becomes 

~j~ + W + SXA^i cos 2t + 8X^^2 cos 4^ + . . .] i:c == 
and 

we have 

It is clear that k is essentially real so long as this series converges. 
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As an immediate application take the equation in Brown s ' Lunar Theory/ p. 107, 

-\-fi^x {l+fm^—- 3%wi^+ (3m^+-y^-w^+--|-m*) cos 2^+^^%^** cos 4^} = 0^ 
at 

where 

PA. 

1—m rz.— n 1-fWi 

then the equation becomes 

j^i'x{l + 2mj,+im,'-^i-^ni'+m^'{3+ cos 2f f^^^mi' cos 4^} = 0, 

which is of the form above^ £ replacing t We may then take 



X' 


m,* 




"^/ 


n'-l. 


64 {2?7ii+ . 


•■)" 


128 



Here m^ is a small quantity and 

\^ 1 + . . ♦ ) 

is of the order m^^, while similarly X^/{n^—iy is of the order rrii^. Also 

^ = ( 1 + mi) { 1 + f m/ ( 1 -^ 2m, + 3m,^) -^%m/}-^ 

= (1 +mO (1 +fwii^— iw^i^+ W-w%^)' 

Thus 

4^/X^ 



^ - 2 



^n { 1 



which is easily seen to agree with the result given by Beown, or by Adams, ' Goll 
Works/ I., p. 187, when we take account of the fact that 

2iKi = 2^*/c in—n^) t ™ 2iK — — nt, 

^ n 

so that, in terms of the quantity denoted by g. 

This example is chiefly useful here as calling attention to the fact that ^^ while 
not exactly equal to 1, is near to it, and consequently the factor x/{ti^-l) is only 
small of the first order in 7%. The same weakness occurs in the terms in ^C"^ ..., in 
the solution. 



DIFFERENTIAL EQUATIONS OF ASTRONOMICAL INTEREST. 137 

§4. In the equation considered by Hill ('Coll. Works/ L, p. 268) the ratio 
4:kiX/{n^—l) is about (2785)~\ and there is a term slightly greater than 4^iX(2785)~'' 
arising in the terms in X**"^^ in the series for /c, in which Ak^X is about 0*5704; and 
the series fails absolutely in cases in whi<ch n is an integer. Then the assumption 
that AC is a power series in X^, and the terms in X which are independent of X, must 
be modified, for reasons above given. The series when n is an integer has been 
considered by Tisserand, ' Bull. Astr./ IX., 1892; modifying his procedure, so as to 
include the case when n is near to 1 as well as that in which n = 1, we may write, in 
accordance with the suggestion of such examples as that above quoted, 

71^ = l + 4X^a + 4X%2+-»-j 
and then, denoting e^^^'^ + e"^**'^ by '^^r? consider the equation 



dt' 



+ [l+4X(/^i + ^ii^i)+4X^(/^2 + /^2'^2) +...]r^ = 0. 



By the changes 

the differential equation may be replaced by the pair 

where 

(p = X{hj + kiWi) -\-X^{h2 + k2W2) + .... 



Assuming here 



q = Xgi + X^g2+..., 
U = 1+Xi^i + X^t^2+..., V = B(l + X'?;i + X^i;3+...X 



in which B is a constant, and u^, Vg are polynomials in ^ and f \ we find, equating 
coefficients of like powers of X, 

_ |_ Q^t(/f._i + ^2^r— 2 I • • • "T Qr — -^r? 

Out 

dv^ 



dr 
in which 



4- git^r-l + g2'?^r-2 + - " + Qr = K 



r> 



H, = {hi + kiWi) (m,_i — f B'y^_i) + (A2 + ^2^2) (Mr-2 - fBv^ ._a) + ... + {h^+ \w^) ( 1 — ^B), 



r 

VOL, CCXYI» — A. U 
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In these equations, as u^^ v^ are to be polynomials in f, ^"\ the absolute terms, 
those involving ^^, must vanish. For r — 1 this gives 

We thus write, using hyperbolic functions. 

With these we find at once by integration the values of Ui, Vj, save for the absolute 
terms in these, which we denote by P^, Qj respectively. The conditions for these are 
to be found by considering the absolute terms in the equations for r = 2 ; and so on 
continually. In general, when we have found 

Ui, Vi^ ^2, 1^25 ..»j '^r-l'> '^r--!? 

and have found u^, v^, save for their absolute terms, P^, Q^, we find, on taking the 
absolute terms in the equations which involve du^^i/dr and dv^^Jdr, and adding and 
subtracting these terms, that the two quantities 

are thereby expressed in terms of known quantities. It is at once seen that there 
would be no loss of generality in putting P^, P^, P3, ... all zero. Carrying out the 
work, and writing M^ for P^— Q^, we obtain 

where - 

hi = kichaj 

h.^ =: M.j^kisha-~-^kich2a^ 

h^ = ^'Mi%e'"-Miki^sh2a + k^^cha {2sh^a-l) + kj,k2Gha + kisha{M2-~''M.iFi), 
Also 

___i^^i(i^2,H ^ l^e-'t+AWi + X'W,+ ... 



in which 



+ P2 — V^sha — M^kicha + kisha {cha + e"-) 

-^e-^ [P2-'M2 + TAsha + Mj€ie^''-h^sha {cha + e""')] 
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If from these formulae we determine Mi and Mg— MiPj in terms of h^ and ^3 we 
find for q, 



q = kiShaX + M.2Ctha . X 



2 



where 



Hs ^ ^ 2JcJl2ch^a h^cha — ki k^ -~ j^/ ( 2sh^a — 1 ) 

2kiSh^a sha sha 



I ' • • • 



This formulae is apparently unsatisfactory when sha is small, or n^—l nearly equal 
to AXki, In fact, the series is of the form 

whose square has a form in which we can put a = 0. On squaring, we have 

q' =: {h,'-h') x'+ 2/hH,x^+x^ {lIi-4hm2+2hA-2h%-K%'+W) + . . ^ , 

wherein 



2 — '^2 'n J 2'n ? 



and this form is appropriate when ot = or hi •= ki. In particular 5 when 
A3 = A3 = , . . = Oj but hi is not zero, this gives 

q' = {hi'-ki')\'+hi{3ki'-^2hi')X'-\-[5{hi'-^ki'^^ 

a formula reproducing the former if h^^-h^ + h^^ be put for Ap It will be seen in 
Part II. of this paper why the form of q^ is comparatively so simple. 

Brief reference may be made to another way in which we may use the foregoing 
equations, regarding h{^ Ag, A3, ... not as given constants but as quantities to be 
determined to simplify the result ; this has been adopted by Prof. W hittakbr 
(' Proc. Math. Soc.,' Edinburgh, XXXII., 1913-14) who chooses as his condition that 
no terms in f°, ^^ shall occur in W^, W^, ..., in the expression for x. This can be 
secured by taking 

Pi = hisha^ Ml = 2kisha, P2 = 0, Q2 = 0, — 

From our present point of view a more natural procedure is to take 

Pj = = Qi = P2 = Q2 = .... Then we obtain 

n^ = l-\-Xkich^-'^X%%h2fi^X^lki^ch^{28h'^-i)+kihch0]+ ... , 

where we have written ^ in place of a, as this argument is now supposed to be 
determined, from this equation, corresponding to a given value of n^. When ^ is so 
determined, q is given by 

q = kiXsh/3-ki\'sh2/3+X'lkM^{6ch'^-^i^) + kA^^ 

V 2 
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an equation which does not contain sh/3 in its denominator. With a view to the 
comparison of this method with the two others given in the present paper we 
consider two examples. First, for the equation 

at '^ 

for which n^ is actually unity, we should determine ^ so that 

= h^ch^-^k^^ch2/3 + Jc^^ch^ {2sh^p-i) + kjc^ch^ + . , . , 

where we. have replaced X by 1. This gives approximately 

ch^ = -ih (1 +W-h), sh/3 - i (1 ~W), 

and hence 

while the value for /3, substituted for a, gives the series for x. We may remark that 
for the equation 

"j^ + (l + 8^1 cos 2^) i^ — 0, 

TissERAND (' Bull. Astr./ IX., 1892, p. 102) finds 

9719 



q — iiCi ll—gki h _^"y ^-g fCi + a . . ) 



As a further example take 

CO JC' 



dt 



+ a? [l + 4^1 (l + Wi) + 4:k%V2+ . . .] = 



which, as will appear, is an interesting equation. Then /3 is to be found from 

^1 =: k,ch/3-^h^ch2/3+k{^ch^{2sh^^-i) +hk2ch/3-^ 



so that 



and hence 



ch^ = 1 + ih + iki - i^s + ih^ - -Pi^s + • » ' . 



In both these examples the value found for q follows at once from the general 
formula above given for g^ of which a further deduction is found below in Part II. 
In the last example the value found for /3 gives a solution for a^ in a series involving 
{ki)K It will be seen in Part II. that when x involves {h)% it is in a very simple 
way, and the case seems better treated as there explained. The occurrence of (k^f in 
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the value of g, in certain cases, is a particular case of Poincare's theorem, * M^th. 
Nouv./ I., §79, p. 219. The phenomenon presents itself, however, as a consequence 
of the use of elliptic functions in Tisserand's theory of the small planets ; see 
TissERAND, ' Mec. Cel.,' IV., p. 426 (or ' Bull. Astr.,' IV.). 

§ 5. A very important question in regard to the differential equation under 
discussion is whether q is real or not, since upon this depends the conventional 
stability of the secondary oscillation determined by the differential equation. We 
have remarked above (§ 3) that when n is not an integer, and ^^X, h^, ... are small 
enough to render the series there obtained convergent, the value of q is necessarily 
real. The cases in which n is an integer and h^^ ~ h^ — ... have been discussed 
by TissERAND, ' Bull. Astr.,' IX., 1892, who obtains the result that the motion is 
unstable for ti = 1 or n = 2, that is for the equations 

-y-g- -f- [1 -f 4XkiW^ X = 0, -jY + [4 + 4x^17^1] a: = 0, 

when X is small enough, but stable for greater integer values of n. The formula for 
5^, given in the earlier part of § 4 preceding, shows that for cases in which 

n^ = l+4AiX 
the motion is stable provided 

ik/hf > 1, 

the values of choL and shoi being then both real. It shows further that it is stable for 

h^ = ±ki = positive 
provided X be small enough. The critical equation is thus 

-\-x[l-{-4:Jci{l-{'Wi)'-\-4:k2W2-\" ..*~\ = 0, 



dt 



TT 



the other sign of ki being obtainable by changing t into t+-^* 

§ 6. We proceed now to the case when n = 2. 
If in the equation 

-yy- + X [m^ + 4X (/h 4- kyW^ + 4X^ ( A^g + k^w^ + . . . ] = 0, 
in which m is an integer, we put 
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we obtain 



where 



dr m dr m 



We may then further substitute 
leading to 



where 






These equations can be solved by writing 






XJi = l+Xt^i + X^t^o+..., W = . -^— +XWi + X'M^a+..., 

m 

where A is a constant, and u^, Uq, ..., w^, w^, ... are polynomials in f, f™\ 

For m = 2, in particular, we find that if h^ = 0, the quantity A is required, and 
determined in the course of the work, and ^i " 0. But if h^ is not zero, we must 
take A = 0, and obtain q^ — ^\, the succeeding q^, ^3, ... being real. In fact, as far 
as X^, 

which gives 

+ { eW + iUi^i - A^i' + Wh - W - K (Pi' + W) + W + PA} ^' + • • • • 

We know, as is shown in Part II. of this paper, that the form of q^ is valid even 
when hi = 0. Then we have 

q' = iX* {h, + Jc,- W) {h-h + W) + ^ * • . 
which, when h^ = 0, is only positive, provided 

i)A/| ^ 0IC2 ^ Kf\ • 
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The case discussed by Tisserand is that in which k^ = Jc^ = ... = 0. Then 

g = ^'A/(-5), 

and the quantity A in the formula for W^ or x^ is found to be 3-( — 2±v — 5). 
When w = 3, for the equation 

-7-^+(9 + 8X^cos 2^)cc = 0, 
we find A = Oj and 

The question of the reality of g, in cases where ^2 = = ^3 = . . . , is discussed by 
PomCARE, ' Meth. Nouv./ 11. (1893), p. 243, and by Callandreau, ' Ann. Observ.,' 
Paris, XXII. (1896), p. 23. So far the results are : — 

(1) For the equation at the bottom of p. 135 (§3) g is real when n^ is not an 
integer, provided the series obtained converges. 

(2) This condition does not however include, for instance, the case when n^ is near 
to unity. For g is imaginary, for the equation 

--TT + [^^ + 8^1 COS 2^ + . . .] ir = 0, 

if {n^—lf < (4^i)^. It is real if {iff —If > (4/^i)Vand real if n^—1 is positive and 
equal to ± 4^i. This has been proved here. 

(3) g may be real when n is just greater than 2, when ^1, ^2, ... are small enough. 
This has been proved here. 

(4) g is real when n is any integer greater than 2, if Jjg = ^3 = , . . == 0, but 
imaginary when n — I or ^ = 2. This result is given by Tisserand and 
Callandreau, as above. "^ 

[October SO, 1915. — ^It may be worth adding, in connexion with the numerical 
results given in § 6, that the equation 

~YT +csm ^. a: = O3 
in which c is small, is solved by 

^ See the note at the conclusion of § 21 (p. 184). 
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in which, as far as c\ 



and, as far as c^ 



A = -7|(1+IV). 



f f sin t + ~~y~ sin 2t - -^--— - 
4/ 31 . . ^-x o^ 11 • Of , cos 4^' 

+ C* 7- ^ COS ^--m COS 2t-— -7= ^ COS 3t + "—r^ 

\8y2 432v/2 4608 

§ 7. We pass now to the consideration of a pair of simultaneous differential 
equations arising in the consideration of the stability of the motion of three particles 
occupying the angular points of an equilateral triangle moving under their mutual 
gravitation. 

The stabihty of this motion has been discussed by Routh (* Proc. Lond. Math. 
Soc.,' VI., 1875; ' Rigid. Dynamics/ II., p. 6l) in the case when the relative paths 
of the particles are circles.^ In what follows we do not assume this. 







The three particles being S, E, M, take an axis through S, say SX, rotating with 
angular velocity 6, the line SE being supposed to coincide very nearly with SX. 
Draw a perpendicular EH from E to SX, denote EH by y, and SH by A -fa;, where 
X, y will be considered small, their squares being neglected, but A is a variable finite 

quantity. Draw a second axes SY through S at a constant angle •- with SX, and 

* The following references may be of use :--Gharlier, 'Die Mechanik des Himmels,' and 'Astr, 
Nachr./ 193, 15; Stockwell, ' Astron. Journ./ 557 (1904); Linders, ' Arkiv for Mat.' (Stockholm), 
IV., No. 20; Brown, * Monthly Notices, E.A.S./ LXXL (1911), pp. 439, 492; Heinrich, ^Astr, 
Nachr./ 194, 12 (December, 1912); Block, ' Arkiy for Mat.,' X., 4 (1914). 
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similarly, draw a perpendicular MK froixi M to SY ; denote SK, KM by A + £ and ?^. 
If Pv = SE, r = EMj p = MS, we have, with proper conventions of sign, 



^^2 __. 



}|(A + ^)~|i;v^3-A-a?J+[|(A+^)\/3+ii|'-|/J. 



The accelerations of E, relatively to S, parallel to SX and parallel to HE, are, 
respectively, 



(E + S) 



A -f" 



R 



X 

3" + 



(A + ^)-|^v/3-(A+a;) _m i(A±£>Z±W^ 



8 



p 



8 



{A+i)y3+h-y ^fh+i{^+i)\/s. 



(E+S) :!, +M^^^ ^^^^\"^^ ''"'^ -M 



R 



3 



^. 



the accelerations of M, relatively to S, parallel to SY and parallel to KM, are, 
respectively. 



(M+S)^f-~E(i.ii^^±^ 



hi \/3 



(A + a?) , v^S- 



1 

2 



( A+^z V 3 + 2"^i— -^ 



r' 



B 



•p |(A+a?)+iyv/3 

li 















E^ 



Ifj then, in the equations of motion relatively to S, after expanding in powers oi 
x^ I/, ^, ??, we equate the finite and the small parts, the squares of x, y, ^, rj being 
neglected, we obtain 



and 



w 



here 



p. 



together with 



M 



2 



A-AO' 



A^(9 = constant, = h, say. 



4" Nj -I- Vi 



A do 
dt 



did „ 
dt^' ' 



X-20Y-eY- 



2 fi 



A» 



X 



O 



o /g 



4 A' 



4S + E4-M 






3^/3 



Y+20X + 0X-f0^~-^JY= ,,, 

AV 4A*' 



^[(E-M)X+y3(E+M)Y: 
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in which X, Y respectively denote ^—x and r}—y^ and also 



x — 26y—Qy— ( 6 



y + 2ex + ex-(0' 



2 M 



w/-JlL. 



/x. 



3 



__ 3v/3M 

X — 



3 



4ya 



1\jL'\/ d 



tc+xys+Y 



\ 



y 



3./3 



4A' 



;x-Yv/3; 



The first equations have integrals expressible by 






= 1 + 2X cos 6, 



If 
I 



fA. 



the point (A, 0) moving in an ellipse of eccentricity 2X and semilatnsrectum /. "With 
these values th^^ other equations are much simplified if we take 6^ instead of the 
time t, as independent variable, as was pointed out to me by Mr. H. M. Garner, of 
St. John's College, Cambridge. With this change they become 



(l + 2Xcos0)(X''-2Y'-X)»-4Xsin0(X'"-Y) == aX+hY, 
(l + 2XcosO)(T^ + 2X^-^Y)--4Xsin0(T + X) = ^X + &Y, 



> • 



^ 



where 



a — 



8S-il 



4/x 



h 



3(E-J 



4/x 



b 



4S+5(E + M) 



ifj. 



and 



de ' 



II 



d\ 



dQ 



2 J 



«c. 



z 



tofifether with 



3, 



{l + 2\(iOBe){x"-2y'~x)-4.\Bine{x'-y)-2x = ^^''^ (Xys+Y), 



•^ 



(1 +2X cos d) {y" + 2x'-y) -^4X sin 6 (y' + x) +y = i?^^^ {X-Y^S). 

4/X 



J 



(I) 



y . (II) 



The first thing then is to solve the equations (I,), after which the right side in (II. 
will be known. Considerable simplification can be introduced by change of notation ; 

1 , /27rA 2 f4:7ri\ 

let w =^ exp[--^j, w ^ exp I ---■), 



A ™ -|(a + &+2) 



3 

2-, 



so that 



XT 



.> S-\-wFj~\~w' 



3 

2" 



H = ^{a—b + 2ih). 



3 B + W^Fj + W 



K^ -\- j2j'{- 



where 



m 



^^ 



S+E + M 






(S-fE + M)' 



2 



{a—b~-2ih)^ 



""x 9/' 
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further 



|) = A+{HK)^ = i 



1+1 






2\f 



q = A- (HK)^ = I 



~9 



2\r 



so that 

Also 



l.v^2 



i) + g=,3, ^g = tml 



tt = (1 4-2X cos 0) (X+iY), t; = (1 + 2X cos 0) (X-iY), 



whereby the equations (I.) become 



(l 4-2X cos Q) {u^^+2m^) = Ate+Ht)^ 
(1 + 2X cos e) {if-2iv') = Ku+Av,^ 



» ■ t » 



• • 



(i-X 



m which u' — du/dd, &c., and then 



so that #5 iP^ are both real, and 

#+f%-2K*(l+2Xcos0(X+^YX *-i%-2H*(l + 2Xcosa)(X~fYX 
and the equations (I.) become 

(l + 2Xcos0)(cr---2%O ==P*,1 

(1 +2X cos e) {%''+2$0 = g^J 

in which, beside the eccentricity 2Xj. there are the two constants p, q, which are 
dependent upon the single number m. 

The equations (II. ), by means of the changes 



become 



U = (l 4-2X cos 6) (x+iy), V = (l + 2X cos 0) (x—iy), 



{1 +2X cos e) {U^^ + 2^'U0 -f (U+ V) 



3M 

2/i 



{l'-w^)v, j 



(l + 2Xcos0)(V'^^2^VO-|(U+V) = ^^ 



2/UL 



■• (l— w)m 



J 



. . / (II y 



Consider now the equations (I.)". We know that the solutions are of the form 

where C, Cj, C2, C3 are arbitrary constants, F, Fj, ..., G^, Gg are definite functions of 
period 27r, and ic, k^, k^, k^ are definite constants. When X = 0, substituting in the 
equations 
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we obtain 

so that the values assumed by k, ki, k.j, /%, when X = 0^ are the roots of the equation, 

or - . 

Thus 

^= ±{i'(l + '?Bp±(l~#'i)''^}, 

and the four values ar^ all imaginary when m > 1 , and all real when m < L 

Supposing S > E > M, we find at once, from the formula for m, that the least 
possible value of S/{S + E4-M) in order that m < 1 is 0'96147..., but thia requires M 
to be very small; but if S/(S + E4-M) be greater than 0'9618..,5 then m is certainly 
< 1 even if E ™ M. In our solar system the sun s mass is more than 99*8 per cent. 
of the mass of the whole system ; thus if S in our problem were the sun, and E, M 
were any two planets of the system, the condition for m < 1 would be easily satisfied. 
We shall then suppose m < 1. 

Now compare with the equations (L)^' the equations 



(1 --2X cos e) {#''-^2^0 == P 



(1-2X cos 6) (-^^^^^ + 2$^ = q^ 



r 



^ .J 



(III.) 



obtained from (I.)^^ by change of the. sign of X. They can also be obtained from (l.y^ 
by changing 6 into 6 + tt. This last change shows that the characteristic constants k 
belonging to the equations (III.) are the same as for {I,y\ while the former change 
shows that the values of k proper to (III.) are obtained by changing the sign of X in 
the constants k appropriate for (I.)^^, When m is such that the values of ^ for X = 0, 
namely, the four values of a- above, are all different, a change in the sign of X cannot 
interchange the values of k among themselves. Thus we infer that each k is unaltered 
by changing the sign of X ; for two of the values of o- can only be equal when m^ = 1. 
In the applications in view of which the question V7as first considered, S denotes the 
sun, E denotes either Jupiter, or another planet such as Mercury, while M is of 
negligible mass. When E is Jupiter we have 

w' = 27r(^5W(l+ToW)'- 0-0257, X - KO'OS) = 0'025, 

and m^lx is nearly unity. When E is mercury 

m^ = 27/5*10^ = 0-0000054, X = 4-(0'2) = 01, 

and m^ = 5'4X^, m — (2'3)X^, nearly. In either case we may regard m as sm.all, and 
the four possible values of cr are approximately 

± (1—- |-m^), ±|-m, 
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of which the j&rst two correspond to a period nearly the same but shghtly greater 

than that of E, and the last two correspond to a period — times that of E. Wlien 

m 

E is Jupiter, this last is ^^ times the period of Jupiter, or nearly 150 years ; when E 
is Mercury, this period is approximately 200 years. As m is small we have approxi- 
mately 

To neglect m^ would be to neglect the ratio 27E/S ; but we may remark in passing 
that if we put q = 0,p = 3, the equations give 

%^ + 2# = C, a constant, 
together with 

1 + 2X cos 

of which the integration can be completed in finite terms. For it may be verified 

that the equation 

(1 4- 2X cos 0) <l>^^+ (l + 8X cos 0) $ =:. 

possesses the two integrals 

sin 0(1+ 2X cos 0), 

cos0-2X(l+ sin^0)--4X^cos + 8X^ cos 20+1 2X^ sin 0(1+ 2X cos 0)yr, 

where 

d6 



^ = \ 



1 + 2X cos 



§ 8. We consider briefly, first of all, what would be the application of the method 
of infinite determinants to the equations (I.)^', which we may now write, with x, y for 
$, %, in the forms 

{l + 2XGOBd){y^^ + 2x^) = qy. 
We should substitute 



GO 



— 00 — 00 

and equate to zero the coefiicients of the various powers of e^^ The substitution 
gives, if f = e^\ 

[l+X(f+^-^)]2:[A^(/c + ^^)H2^'(/c+?^)BJ^+p2A,,^ == 0, 
[l+X{i+t~')']^[B„iK+nf-^2i{K+n) AJ ^+g2B,r - 0, 

and, denoting K + n by k^,, we obtain for the unknown coefficients A,,, B,, the equations 

X ( A^_ia:\_i + 2?.'B,,_i/c^„i) + A^ (/c/ -i-p) + 2^"B,,/c,, + X ( A,,.^i/c\+i + 2iB^,^jK,,^j,) = 0, 
^ (■--2^A^_iK-^„i + B^„i/cVi)--~2^*A,,/c^ + B,,(/r/ + g) + X (-2^A^^ ^ 
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If we now write 



P = 



^Kn +2^B„/C 



n'^n? 



^?n 



2?/A„/C,, + -0^/C„ 



n'^n ' -^-^n'^^n ? 



which are equivalent with 



A 



K^¥.. — 2i 



^n-^ n 



n 



,'^n('^/-4) 



2iF^ + /c^Q 



n 



n 



^n \^n ^ / 



the equations may be replaced by 






> 



. . (A) 



wherein 



a. 



n 



1 + 



P 



2^p 



K. 



n 



n 



^n \^n '"*/ 



C„ = 



n 



2^'g 



dn^ 1 + 



g 



'^w 



-4 



so that 



fn*-f„' + im2 






it being remembered that p + Q^ — 3, jpg = Jm^. 

When we eliminate P„_i5 ...5 Q„^i from the equations (A), we obtain an infinite 
determinant, which, leaving aside questions of convergence, we may denote by 



a 



-1 



X 



A c_i (i„i 



X 



X 



X 



a 



c 



X 



X 



d 



X 



ai 



X 



X Ci d^ 



X 



0. 



The product of the diagonal determinants cbn'^n'^^^rfin i^ ^^^^ 

sin 7r(K:— 0-1) . sin tt (/c— o-g) . sin tt {ic—a-^ , sin tt (ac—o-^) 



sin^ TT/c 



where o-i, erg? 0-3, 0-4 are the four roots of (T^—cr^-\-^7Y}? ™ 0, previously considered. In 
using this determinant to obtain a further approximation to k it would seem 
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appropriate to use a theorem* for the expression of a determinant of 2n rows and 
columns as a Pfaffian, a sum 1 . 3 . 5 ... (2n— l) terms, of which each term is a product 
of ti factors, each factor being of the form 



(12) = aib\-a\b^ + a2b^ 2- (^^2^2-^ ••• + <^J^^ n- <^^nK 



where the elements 






/ JJ 



€v jO J J (a/ 2^ 2j * * * Co ^0 



n> 



are the constituents of two rows of the determinant. For in this case the factors (12) 
are easily calculated. But we do not pursue this method. 
§ 9. Instead we proceed as follows. In the equations 



where ^ = e^^, write 






in which Kg, k^, ... are certain functions of p, q to be determined. Then the equations 
become 

[i+\(^+ri)] [x"-2Y'+2*k(x'-y)-/c=x] =px, 

Cl+x(f+r^)] [Y"-2X'+2*.(Y'+X)-.^Y] =. qY, 

which by the general theory are capable of periodic solution when k is properly 
chosen. Put then 

.X. = 1 + X^>i + X^2+...5 x=.r(l + Xi^i + X ^f/'2 "^" • • • )' 

where P is a constant ; the differential equations then take the forms 

{l + \w) {H.Q + \tii + X ti2+ ...) =jpAj 

(1+Xw; (Ko + XKi + X K2+ ...) = gi, 
w denoting ^+^~^ Comparing the coefficients of like powers of X^ 



and, in general, 
so that 






* Proved in Sgott-Mathews' 'Determinants' (1904), Chap. VIII, p. 99, § 19. Also in Baker, 
' Multiply-periodic Functions,' p. 314. 
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and, in general, 

where H^, K„ are the coefficients of X^ respectively in 

X^^ ^ 2 Y^ + 2^'/<: (X^ - Y) ™ k'X, 
T' + 2X^ + 2u ( Y^ + X) ■»- k'Y, 

Ho =- -2^VP--(T^ Ko = 2^0" -o-T, 



In particular 
so that 



and, as previously, 
while, if we write 



o-^+J9 + 2^VP == 0, 2^cr ==:: [a-^^q)V, 



— 2^cr 2^(r 



which are both pure imaginaries, we have PQ =: L 

Next 

Ki - P [VA + 2Q/i + 2*0- (VA + Clin) -<7=V.J ; 

putting these respectively equal to p{<j)^ — iv), Vq{\j/-^—w), we obtain two differentia] 
equations for (p^ and yjr^. If we assume 

and notice that 

we find, writing o-^ for a- + ?^, 

Ai (a-/+_p) +2P^ViBi ^ p, A„i (a-„/+p) ~f 2P^V„iB-i ^ P. 
—Ai . 2Q^Vi+ (o-/ + g) Bi =^ g, — A_i2Qicr„i+ (o-_/ + g) B„i :-^:: g. 



If 



4 2 1 1 2 



these give 



-(o-/-fi3)g+-H^'+g' 



with similar equations for A_.^, B_3, 

Proceeding similarly to equate terms in X", w^e find 
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If herein we assume 

and equate terms in f^, ^""^, ^^ we obtain 

A2 {0-2 +p) + 2P^'cr2B2' = p (Ai- 1), A_2 (cr^_„2+.p) + 2P^V__2B_2 = ^ (A_j - 1), 

"A2.2Q^V2+((T/+g)B2 = g(Bi-l), -A„2 . 2Q^V„2+ ((t^2+V)B_2 - g(B_i-l). 

and 

(cr^+jp){H — K)+2/Ci>(^P + (r) =_p{Ai + A^i-2X 

-^{cr^ + g)(H^K)+2/c2(-^'Q + (T) = g(B, + B^i-~2X 
wherein the coefficients of H — K and K2 have for determinant 

which is 

and is not zero. That H, K should not be determinable separately is obvious 
a priori ; to regard H as zero would be equivalent to dividing X, Y by a power 
series in X^ with constant coefficients. We notice that the successive coefficients 
Aj, A^i, ..., B2, B„2 are all real The value found for /cg is 

/c2-pg 



cr{l-2u^'){l-4.a-'Y 



A similar procedure can be continued. The differential equations for ^g, -^3 can be 
solved by forms 

the differential equations for ^4, i/r^ by forms 

^, = A,r+A_,^HM2r+M^,r^+M, 

and then the terms in f ^ will involve the unknown quantities 

(o-^-fp)(M-^N)-f2/c4(^P + 4 

-~(a-^ + g)(M-N) +2/c,{-^Q4-crX 

from which /C4 is found. And it serves as verification of the computation to see that k^^ 
involves H, K only in the combination H— K, as it must in order to be determined 
without ambiguity. 

VOL. CGXVI.' — -A, y 



154 PROF. H. F. BAKER ON CERTAIN LINEAR 

The value found for x is of the form 

or, say. 



X = e 



where every one of ui, u^, u_^, u^, u__^, ... is a power series in X^ with real coefficients, 
not generally vanishing with X\ And similarly for y, 

§ 10. The interesting case of the preceding solution is that corresponding to the 
value of (T given by 



The quantity 



is then equal to 



cr = |-[{l+m)*--{l™m)*l =im{l + '^ + ...). 






^m (7 + ^mP) 



approximately, and /c^X^ is of the order mX^ When ?7^^oc x, this is of the order m^ or 
^hi2 . ^^]^Qj^^ r^-^ oc x^ it is of the order m^^'^ or X^ Thus a very few terms of the preceding 
solutions would seem to be sufficient for practical cases. 



PAiiT ir. 

§11. A large part of the interest of Poingare's ' Methodes Nouvelles de la 
Mecanique Celeste' depends on his criticism, of the convergence of the series used 
by astronomers, particularly those series in which the time enters only under 
trigonometrical signs. In t. IT, p. 277, he refers to a linear differential equation 

■s 

^+i»(l+^) = 0, 

in which ^, for our purposes, may be supposed to have a form 

yjy^ c™ 4a cos ht + ih cos M^ 

in which a, h are small. When h^ h are commensurable the equation has periodic 
coefficients, and Poingare makes the convergence of the series expressing the solution 
depend on this circumstance {' Meth. Nouv.,' t. I., p. 66), Considering the case in 
which h and k are incommensurable, and so \// not periodic, and supposing a, 6 to 
have common a small factor /x, he obtains formal solutions of the differential equation 
in sines and cosines^ and says ^' les series . . ., qu'on pent ordonner suivant les 
puissances de fi, ne sent plus convergentes " ('Meth Nouv.,' t. II., pp. 277, 278). On 
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the contrary, I believe that the solution of the differential equation above, arranged 
as a power series in a and 5, converges for all finite values of these parameters, 
and that this is a consequence of a general theory of linear differential equations 
considered in papers^ , published by me in 1902. As this theory is capable of 
application to many other differential equations, as will be illustrated below by 
application to the equation considered by G. W. Hill for the motion of the moon's 
perigee, I wish to deal with it here, repeating the argument in part. 

§ 12. Consider any system of linear differential equations, the v? coefficients 

/# ''V* 

~ « — ^^il 1 ■" • • • ^1 ^XYi^Yi^ \ ^ J-J ^5 •••^ /? 

at 

u^j being functions of t. If these are considered only for real values of t, the 
properties which we require to assume are that, along a certain range which we 
shall suppose to include t ~ 0, these functions Uij are single- valued, limited, and 
capable of integration, the same being true of certain other functions derived from 
these by multiplications, and further, that certain infinite series, which we shall 
prove to be absolutely and uniformly convergent, are capable of differentiation, term 
by term. But in the majority of practical cases the coefficients u^j may be looked 
upon as the values, when t is real, of functions of a complex variable t. In this case 
we suppose a star region to be defined by lines passing to infinity from certain points 
in the finite part of the plane, which we call the singular points ; we suppose t ~ 
not to be a singular point, and the lines may be straight continuations of the radii 
joining the origin to these singular points. Within this star region, bounded by the 
lines in question, the functions u,j are supposed to be single- valued and capable of 
development by power series about every point, forming monogenic analytic functions 
in the usual sense. Taking then any region within this star region, we obtain 
solutions of the differential equations, with arbitrary values for ^ ~ 0, in the form of 
infinite series of functions, obtained by quadratures, which are proved to converge 
absolutely and uniformly within the region taken. 

The method of forming these solutions is extremely simple, involving only 
integrations and multiplications, but the way in which the work is arranged, though 
often of great utility, does not seem yet to find common acceptance, and some words 
must be given to it. 

* 'Proc. Lond. Math. Soc.,' XXXIV., 1902, p. 355; XXXY., 1902, p. 339. See also the same 'Proc.,' 
2nd Series, IL, p. 293, where it is explained that the same idea had already been used by Peano 
and others. To me the method was independently suggested by the theory of continuous groups, 
'Proc. Lond. Math. Soc.,' XXXIY., 1902, p. 91. Poingare's conclusions as to the convergence of 
astronomical series have been criticised by G. W. Hill, ' Coll. Works,' IV., p. 94; but the point there at 
issue is different from that considered here. In connexion with an example considered by PoiNCARli:, 
loc. cit., p. 279, see Bruns, ' Astr. Nachr.,' No. 2606 (GIX., 1884), pp. 217, 218. Also Borel, ' Theorie 
des Fonctions^ (1898), p. 27; Hardy, ' Quarc. Journ.,' XXXVI., p. 93: ^ Proc. Lond. Math. Soc.,' III., 
p. 441, and the references there given. 

Y 2 
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The v? quantities Ui-j can be arranged to form a square of n rows and n columns, 
the first suffix % denoting the row, and the second suffix y denoting the column in 
which a particular element u^^ is placed. This square is denoted by a single symbol, 
say u^ and called a matrix. The symbol %iv, formed from the two symbols te, i;, 
written in a defini|;e order, denotes then another matrix whose (^, y)*^ element has 
the value 

which is formed from the elements of the i^^ row of the matrix u and those of the 
j^^ column of the matrix v. This new matrix uv is called the product of u and v, 
taken in this order ; it is generally different from vu. The symbol 1, when used for 
a matrix of an assigned number of rows and columns, denotes the matrix of which 
every element is zero except those in the diagonal, all of which have the same value, 
unity ; it is easy to see that any matrix is unaltered by multiplication with the 
matrix unity of the same number of rows and columns. The symbol w'^ denotes 
the matrix such that the product u~^u is the matrix unity ; in that case uu~^ is equal 
to ir^u ; the symbol u~^ is nugatory when the determinant formed with the elements 
of u is zero, and only then. In general, the determinant formed with the elements 
of u will be denoted by \u\. By the 5t^m, u-\-v^ of tw^o matrices u^ v, of the same 
number of rows and columns, is meant the matrix whose {i, JY^ element m u^j + Vij, 
and, similarly, for the difference. Frequently we denote the aggregate of a row of 
n quantities, x^^ x^^ .-? ^n by the single letter x ; then if ii be a matrix of n rows and 
columns, the symbol ux denotes a set of n quantities of which the ^^^ is 

By the differential coefficient of a matrix we mean the single matrix whose elements 
are the differential coefficients of the given one. In what follows, if the (^, jY^ 
element of a matrix '^ be a function of t, we denote by (^u the matrix of which the 
{i^jY^ element is the integral of u^j taken in regard to t from ^ — to i^ ^ ^. If, for 
an instant this matrix Qu be denoted by v, the product matrix uv will be denoted 
by uQu, and the matrix Q {uv)^ or Q (uQu), will be denoted by QuQu. Similarly, 
Q {u . QuQu) will be denoted by QuQuQu, and so on. 

Now consider a matrix of which the (i, jY^' element is the infinite series formed by 
the sum of the {i, jY^^ elements taken from the matrix unity (of the same number of 
rows and columns as u), the matrix Qu, the matrix QuQu, the matrix QuQuQu, and 
so on. This will be denoted by 

Q (u) ~ 1 + Qu-\- QuQu + QuQuQu + . . . , 

and the series on the right will be said to be uniformly and absolutely convergent 
when this property is proved to hold for each of the n^ infinite series which constitute 
its elements. 
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Repeating now the demonstration given, ' Proc. Lond. Math. Soc.,' April 10, 1902, 
p. 354, let 7ij/^' denote the {i,jy^ element of the matrix Qu, that is. 



«,« = 



W^'j U/L y 



'IJ 



similarly, let u^j^^^ denote the (^, jY'^ element of the matrix QuQu, namely. 



Uij 



(2) 



H 




lU^lUy -jr y^i2^2j + . . . + y^in'^nj J ^^J 



and so on. For the region chosen within the star region above explained, when the 
fmictions u^j are functions of a complex variable, or for the range of values of t 
adopted when the elements u^j are functions of a real variable, there v^ill exist a real 
positive quantity M^^ not exceeded by the absolute value of u^j for the values of t 
involved. Taking a path of integration limited to such values, from the origin i^ = 
to t = t, this being a rectifiable curve of length 5, let ti be an intermediate point of 
this path, the length of the path from the origin to ti being Si, Then we have, 
considering absolute values, 



and in particular 
Similarly, 



Jo 



"^ij \h) I -< "^i-^Vl^--. 



tj' 



u/^ (t) I IS (M,i5iMi^. + . . . + M,,sMnj) ds, ; 



now if M denote the matrix whose (i^JY^^ element is M,;^, the (^, ^'f^ element of the 
matrix M^, formed by the product of M with itself, will be 



-il^y-i-lj 



.,-2XTj-2j 



m-^'-^nj > 



which we may denote by (M^),;^ ; thence 



and in particular 



u,mt)\ = (M.% s,ds,=is^M%, 



uf (t,) I = W {M.%. 



We can continue this process. The next step will be 



'^ii^ (015 W ds, [M,, {M.%j +...+ M,„ (M^),,], 



3! 



mij. 
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Thus we see that each of the n^ infinite series constituting the elements of the 

iXlOJ UX X.A, 

Q{u) = l-{'Qu-}-QiiQu+ .^^ 

has terms whose moduli are respectively equal to, or less than, the real positive terms 
of the corresponding infinite series constituting the elements of the matrix 

2! 3! . 

This last is, however, certainly' convergent for all finite values of s, whatever be 
the (finite) values of the elements of the matrix M. For the case when the algebraic 
equation satisfied by M has unequal roots, its sum is given by the formula of ' Proc, 
Lond. Math. Soc.,' XXXIV., February 14, 1901, p. 114, which can be easily modified 
to meet the case of unequal roots. 

Thus each of the elements of the matrix Q{u) is an absolutely and uniformly 
convergent series; in the case when the elements ti;^^ are functions of the complex 
variable, as explained above, it follows that every element of the matrix (u) is a 
function of the complex variable, and differentiation (and integration) of the series 
representing this element is permissible, term by term. For the case of real functions 
we introduce this as a condition. 

Hence, if i^t^V denote a row of n arbitrary values x-^\ a?/, ..., x/, ^^^e row of n 

quantities denoted by 

X = Q {u} X 

is at once seen to form a set of n integrals of the differential equations, reducing for 
^ = to the arbitrary values x^, that is, x^ reducing to x^. For if v denote any 
matrix, of n rows and columns, whose elements are differentiable functions of t, if x^'^ 
denote a row of n constants, and y the set of n functions given by 



that is, 
we have 



y z=z q^)X^^ 



Iji — '^•il^l » '^'i2'^'2 1 • » • r ^inr^n } 



dyi ___ dv,^ ^ , . ^^^ 



in ryt t> 



a^ at at 



n > 



which, if ^ denote the matrix whose elements are ,the differential coefilcients of the 
dt 

elements of t% we can denote by . 

dy dv A 
at dt 

Hence the equation 

/ \ '/v 

Ji/ <^ — 4^ \ ^ / *^' ? 
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gives 



dx 


d 


dt 


dt 




_ d 



dt 



\ + Qt6 4- Q (uQu) + Q (uQuQu) + . . . 



X' 



[u-\-uQu-\-uQuQu+ ... 
u [^l + Qu+QuQu+ ... 






x 







or 



dx/dt = ux, 



so that the functions x = Q (u) x^ satisfy the differential equations. By the definition, 
QUij reduces to zero for ^ = ; hence {u) reduces to its first term, the matrix unity, 
when ^ = ; that is, cc = (u) x^ reduces to x = x^ when i^ = 0. 

In what follows we shall require a particular property of the matrix O (u), given in 
^ Proc. Lond. Math. Soc.,' XXXV., December 11, 1902, p. 339. If % v be any two 
matrices of n rows and columns, of similar character to the u considered above, the 
property is expressed by 

Q (u + v) = Q (u) Q {[Q (u)]--' vQ {u) 



where [Qiu)^"^ is the matrix inverse to Q{ii), defined above, such that f 0(?/)]~-^ ii{u) — 1. 
The theorem is nugatory when the determinant of Q (u) is zero. It is only equivalent 
to eayine^ that if in the system of linear differential equations 



that is, 



dx 
'dt 



{u + v)x, 



(jjX^ 

~dt 



(Uii + Va) x^+ .,.+ {Ui^ + Vi^) x^, 



we introduce a set of n new dependent variables, denoted by z^ by means of the 

equations 

X ~ Q (u) z, or z ^ [0 {u)Y^ x^ 

then 

dzjdt ^\Q{u)Y'^ vQ{u)z. 

This follows at once from 



{u-\~v) 



Ju — 



dx ___ d 
dt dt 



[Q {u) z] ~ ~Y Q {u) z-i-Q (u) 

dt 



dz 
dt 



uQ iu)]z-\-Q (u) -r- = uQ (u)z-^Q (u) -r- 
■ ^ ^-^ ^ ^ dt ^ ^ ^ ^ dt 



= ux~{-Q{u) 



dz 
dt 
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which gives 



^ ^ dt ^ ^ 



In what follows we shall generally write 0"^ {u) in place of [0 ('^^)]~^ 
Another property to be noticed"^ is that the determinant of the matrix {u) is 
equal to the exponential of the sum of the integrals from to ^ of the diagonal 
elements of the matrix u. For, if O^j denote the general element of Q{u), the 
equation 

already remarked, is the aggregate of the equations 

Further, the differential coefficient of a determinant of n rows and columns can be 
written as a sum of n determinants, each of which is obtained from the original 
determinant by replacing the elements of one row respectively by their differential 
coefficients. Hence we at once see that, if A denote the determinant of Q (u), 

which establishes the result in question. 

In particular, if the sum of the diagonal elements of w, 



'^11 • ^22 > • • • ~tU 



nm 



be zero, the determinant of Q {u) is independent of t, and is thus equal to unity. 
This result is of frequent application. 

§ 13, After these introductory remarks we may at once show that the equation 

-~^ + xil + 4:a COS ht + 4:h GOB U) = 0, 
dr 

to which reference has been made, is capable of solution as an absolutely and 
uniformly converging series in a, 6, whatever h and k may be. It will be as simple, 
and of utility for other applications we wish to make, to take the equation 

in which we may suppose n to be an integer. 

* Gf. Darboux, 'Compt. Rend.,' XO. (1880), p. 526. 
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111 this last equation, put 



leading to 
then we have 



Writing 
these are 



X ~ 1-6"'^ I -z — inx j , Y ~ |-6~'''^ I -J- + inx 

-inx ^ Xe-'""^ -Y e'''\ dxjdt =- X6-^'^^ + Ye"^^ ; 

dt 2n ^ ^' 

dt 2n ^ ^' 



± (X~Yr), ^ - -^ ^ (Xr-^-^Y). 



dX^ „ 

dr 472. ^ " ^"^ ^' c?T 4.71 



or, say, 



where^ as is usual, the single quantity — •^, written before the matrix, is to be 

multiplied into every element of the matrix. 

In particular, when n =^ 1^ \fr ~ ia cos ht + ib cos kt, 

f-(X,Y)=:(«p + &g)(X,Y), 
where p, q denote the matrices 

-r\ 1/ \-r\ 1 

Thus the solution is expressed by 

{X,Y) = Q{ap+bq){X\Y% 

where {ap + 6g) is of the form 

1 + aQp + hQq + a^QpQp + ah {QpQ,q + QqQp) + b^QqQq + . . . , 

and we have proved that this series is uniformly and absolutely convergent. 

If we assume such a form of solution it is easy by successive steps to obtain the 
values of the coefficients independently of the method we liave adopted. What is of 
present importance is that we have shown the series to be convergent, a fact which 
appears to be denied by Poinca.iie. 

VOL. OOXVI.^ — ^A. Z 
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§ 14. Leaving aside this point, we pass on now to the apphcation of the general 
method here explained to the computation of the inteOTals of particular differential 
equations with periodic coefficients, as, for instance, the equation for the motion of 
the lunar perigee, considered by G. W. Hill, 

It is known from the general theory that the solutions of the n equations 

in which %i, ..., t% are single-valued functions with a common period, say w, can be 
written, in the most general case, in the forms 

wherein Aj, ... , A,, are arbitrary constants, Xi, ... , X,, are n definite constants, and the 
functions (j>ij are n^ definite functions all with the period tv. In many applications it 
is the constants Xi, ..., X^ which it is of most importance to find; when these are all 
pure imaginaries, the motion'^' represented by the differential equations presents, 
beyond the fundamental period w, secondary oscillations of periods 2iTli\, and the 
motion is conventionally said to be stable. 

We show first how this form of solution naturally arises from the point of view we 

have adopted. 

Write Oo* {u) in place of Q (u), and for simplicity write only two rows and columns 
of the matrix, though the argument is quite general. Make the limitation, which, 
as is well known, does not cover all cases, that there exists a matrix of constants, 4, 
of n rows and columns, whose inverse is denoted by /^""^ such that the complete 
matrix Qq"^ {ii) can be written in the form 



Q^'-{u) = h/e^'^'% \h 



-I 



0, e''^"^ 



with only diagonal elements, here denoted by e''^'^ e'"'"^, in the reduced matrix. This 
will be so, in the technical phraseology, if the matrix Qo'''" {^t) has linear invariant 
factors. t Then, from the definition of Q{u), 

while, as u has period w, 

* Interesting physical examples are given by Lord Rayleigii, ' Collected Works,' III., p. 1. 

t A proof of the general theorem for the reduction of a matrix, valid when this is of vanishing 
determinant, is given, ^ Proc. Camb. Phii. Soc.,' XII. (1903), p. 65. The literature of this matter, which 
begins with Sylvester, 'Coll. Papers,' I., pp. 119, 139, 219, and Weierstrass, <Ges. Werke,' I., p. 233, 
is very wide. The reader may consult Muth, * Elementartheiler,' Leipzig, 1899. 
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H 



ence 



and so 




Qo""^' {u) . A = Oo* ('^) . h/e''^"^, \h-\ 

This shows that the matrix on the right has period w. Put then 

which has period w, and is such that V^ = P^^ = 1. The matrix Q^ {u) can therefore 
be written in the form 

0, G'''\ 

which is the theorem in question. 

We now compa^re this with the form of sohition of the original differential 
equations by the method of successive approximation followed by Lagrange, 
Laplace, and others. We have 

\ / ^ 2! ' ^ 



0, e^'^V \0, ^2/ \ 0, c^ 



thus 



1},^ [u) - Po^ + ^P {hyh-^) + ~ P {hy'h-') + . . . , 

Zi 1 



where P is written for Pq*, and y is written for 

Vci, 0\. 

If then, as in Laplace, ' Mec. Cel.,' Liv. II., Ch. V., t. I., of the edition of 187-8. 
p. 266, we obtain the solutions of the differential equations in the form 

where A, B are certain periodic matrices, and x^ is a row of arbitrary constants, 
we can obtain the constants ic^, ic2, tvhich are the most important quantities in many 
applications, hy taking the matrix A, which arises as the coefficient of t^^' and is equal 
in our notation to 

z 2 
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putting therein t ^ 0, so obtaining^ say Ao, equal in our notation to hyh"^, and then 
solving the determinantal equation 

Ao--X| ^~ 0, 

ivhose roots are ie^ and ic2. This process will be found to be equivalent to the 

general procedure explained by Laplace, in the passage above referred to, for 

bringing the time under trigonometrical signs. We have considered only the case 

of linear differential equations with periodic coefficients, and have supposed Qq"^ (^0 "to 

have linear invariant factors ; Laplace's method, if less definite, is of much wider 

application. An interes-ting exposition of the method in general is given by 

M. O. Gallandreau, 'Ann. de TObserv. de Paris,' XXIL, 1896, pp. 16, 20. 

We may notice that 

Aq ™ h/ie^, \h~''^ 



gives 



so that we also have 



0, ic2, 



fi(Ao) = h/e''^\ \h-'', 



0, e''' , 



Q (u) = Po' . li (Ao) 

=: Pq (l+^Ao+2'^ Aq +-"«)3 



and the quantities e'""'"^, e'''''^ are the roots of the equation 

Qo'^{u)-p\^ 0. 

§ 15. When the sum of the diagonal elements of the matrix u is zero, the 
determinant of Q (u) is unity, as above remarked. In this case, when n ~ 2, the two 
quantities e'''''^ e''^'''' are inverses and c^ == — c^. In this case the equation 

Qo"{it)~-p\-=^ 

gives at once the value of cos ew. This appears, however, a less advantageous way 
of determining c^, Cg than that explained above, as requiring greater approximation 
in the calculation of Oq^ {u), as will be seen in examples. 

The fact that Cj, c^ are equal and of opposite signs is a particular case of a 
well-known theorem for the variational equations arising in the general dynamical 
case, which is proved by PoiNCARE (' Meth. Nouv.,' I., 193). The following proof, 
though longer, appears more fundamental in character. The general dynamical 

equations beihg 

dx^ _ _3F dy^ _ _ _9F 

dt dy^ dt dx^ 
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where it will' be sufficient to suppose r to have the values 1, 2 ; let 

be a solution of these equations. Substitute in the differential equations 
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x^ = cj>, {t) + i,, y, =- Va- (0 + nr^ 



and retain only first powers of the quantities ^^ and ri^, which are supposed to be 
small. We thence obtain a system of linear differential equations of the form 

where /3 is the skew -symmetrical matrix of constants given by 



^ 







1 







































1 







(so that /3 ^ = — /3), and A is a symmetrical matrix whose elements are functions of t. 
We then have the theorems following : — ■ 

(a) The roots of the determinantal equation for X, 

|/3-iA-X|=0, 

fall into pairs of equal roots of opposite sign ; 

(b) The determinantal equation for p, 

is a reciprocal equation, unaltered by changing p into p~\ 

To express the proof we require a notation for the matrix obtained from a given 
matrix u by interchanging its rows with its columns, thus placing the element u^j in 
the (y, if^^ instead of the {% jj'^' place. This transposed matrix may be denoted by 
trs (ti) or by u. It is easy also to show that 

[Q(tO]~' = trs[Q(-^)]. 

Then (a) is immediate from the obvious relations among determinants expressed by 



A-~-/3X|-= |A-~^X|-: IA + /3A 



since A — A, ^ = ~^. 
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For (6), since 



/3Q (u) /3~^ - Q {^u^-'l trs (^-^A) = ^A^-\ 



we have the following transformations of matrices 

]Q (/3-^A)]~i - trs [Q ( A^"^)] - trs [f3Q (/3-^A) ^-^] = ^~' [trs Q (^"'A)] ft 
and hencej writing Q~^ (u) for [Q (^0]~^ "th^ following equations among determinants 

0-^ (,3-^A) -p I - I trs {^"^A)-p I = I Q (/3-^A) -~p I , 

which establishes the result in question. 

§ 16. In many dynamxical applications the matrix A is a sum of two matrices 



A 



a 



+ S-, 



where a is a symmetrical matrix of real constants, and & a symmetrical matrix whose 
elements are small. Suppose, further, that p, denoting a row of 2n real variables 
Pu P2y • • • ? "the matrix a is such that the quadratic form 

does not vanish unless every one of the 2n elements of p is zero, which requires that 
the determinant |a| is not zero. Then, if this quadratic form be denoted by ap^, 
and if each of ^ and ^y be a row of 2^^, real quantities, the form 



a{^+if]) ii—in), == a^^ + ia{r]^—^Tj) + atf, 



( 



C^ M iA^\ 



ay^ -f-rj 



has the same property. 

When this is so, it can be shown that the roots of the determinantal equation in ^y, 

are pure imaginaries, and that the invariant factors of the matrix ji'^a—yfr are linear. 
As the proof is not long it may be given here {of, ' Proc. Lond. Math. Soc.,^ XXXV., 
December 11, 1902, p. 380). 

Let -^ satisfy the determinantal equation 

a — /3i/^| — ; 

as the determinant |a| is nol^ zero, \lr cannot be zero. Then 2n quantities x^^ x^, ..., 
whose aggregate is denoted by x, can be taken to satisfy the 2n linear equations 

{a — ^\Jy)x == 0. 

If Xq denote the row formed by the 2n quantities which are the conjugate complexes 
of those of X, we have in turn 

oLXX^ ™ \fr^xxQ, oIxqX ■=■ 'xJyfSx^x.^ ax^x — —yfj^^x^x, 
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and, therefore, -v/^o being the conjugate complex of \fr 

Hence 

(- + —) axx^ ^ 0, 

showing that xl^ + xf^^^ ™ 0, which proves that -^^ is a pure imaginary. 
Writing X for \fr~^, the equations above are the same as 

(a-'/3-X)x = 0; 

we prove that the invariant factors are hnear by showing"^ that it is not possible 
to find a row of 2n quantities y/i, y2^ . . . , such that 

{oL~^^—\)y — X. 
For this would involve 

{^—a\)yxQ = olxXq^ 

of which the right side is real, so that, X being a pure imaginary, either of these 
would be equal to 

(/3 + aX)^o^; = {l3 + Ei\)xyQ = (— ^ + aX) q?|/o, 
of which the last is zero in virtue of 

As axxQ is not zero, the assumed equation for y is impossible, and the invariant 
factors are linear. 

From this fact it follows that it is possible to find a matrix h such that 

h-'^-'ah = / ^0-1 

-*<Ti 

io-2 

-i(T2 

where a-^, a-^ are real. Then the given differential equations, which are of the form 

if transformed by the linear substitution 

iii^ n\y ^2) ^2) = h (Xi, Yi, X2, Y2) 

^ See, for example, ' Proc. Camb. Phil. Soc.,' XII. (1903), p. 65. 
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take the forms 



rfXi dii 0X3 di2\ _ / , o\ /Y V Y X \ 



where a- denotes the matrix above written, with only diagonal elements ia-i^ &c. 
and 9 is the matrix 

The solutions of these equations are then expressed by 

(x„ Y„ x„ y,) = G(<r+e) (Xi», YA x/, y/), 

where X/, Yi°, ..., are the initial values. Now, by a previously given formula, 
where Q (a-) has the simple form 

Q (a) = 



^ia-it 








\ 





g-io-,« 














pio-'it 














g-ior.« y 



the solution is thereby expressed in powers of the small quantities occurring in S-. 

The preceding work has wide applications ; a particular case is that of the 
oscillations of a dynamical system about a state of steady motion, for which S-, and 
9, is zero. 

October 30, 1915. — ^To prevent misunderstanding, two remarks may be added to 
§16. The condition that the quadratic form ax^ should be positive, though sufficient, 
is not necessary in order that the roots of the determinantal equation (/3"^a— ^r) = 
should be pure imaginaries. For instance, if a, &, u, v be real positive constants, and 
H be a quadratic form 

H = Ja {pi-nx.,)' -\~^h {y^-mx^Y - ^ .^r/- 1^ xi, 

the motion about x^ ^ 0, .Tg == 0, ^1 = 0, ^2 ™ expressed by the equations 

ii rr dH/dy^, y^ ^ -dH/dx^, x^ ^ cRjdy^, y^ = -3H/a^2, 

is instantaneously stable if ah{m-nf> (^^ + v)^ the corresponding quartic equation 
having all its roots purely imaginary. This essentially is the case noticed by 
Thomson and Tait, ' Natural Philosophy/ 1., pp. 395, 398, where the illustration is 
that of a gyrostat balanced on gimbals. A simple illustration is also that of the 
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oscillations about steady motion of a weight suspended by a string of which the 
other end is made to describe uniformly a horizontal circle, in the case in which the 
string intersects the vertical drawn dowmvards from the centre of the circle described 
by its upper end. This motion is not, however, secularly stable when there is 
Dissipativity (Thomson and Tait, as above, p, 388) ; and, of course, not instantaneously 
stable, the roots of the corresponding quartic equation having real parts of which some 
are positive. 

A second remark relates to the generality of the form of the differential equations 
used in the text. Equations such as 

dt\dxj cx^ cx^ dx,, 

where ^^^ is a function of Xi^ ..., x^^ capable of expression in terms of n functions 
/3i, ..., ^^ in the form. 

GXg Xy, 

are included in this form, with a slight modification due to the presence of the 
Dissipativity F, and the supposed non-conservative forces Q,,. For this it is only 
necessary to take 

. dJ-j • oJL -T- 

and to eliminate x^, ..., a;„, in the familiar way, from the n equations 






Then the final equations are 



. m . _ _sH_aF ^ 

cy^ dx^ oXj, 

Particular illustrations are : (l) the equations of Thomson and Tait {as above), p. 392, 
for which the coefficients /3,,^ are constants. Then we may take ^^ == c,.iXi + ... +(3,.,,^^, 
where the constant coefficients c,,, are in part arbitrary ; (2) the equations of Lord 
Kelvin for liquid motions of ring-shaped solids, ' Collected Papers,' IV. (1910), p. 106 ; 
(3) the equations of motion of a system relatively to a rotating frame (Lamb, ' Hydro- 
dynamics,' third edition (1906), p. 294. Cf. Thomson and Tait, as above, § 319, p. 307, 
and p. 319), for which we may take, if (f, rj, f) be the co-ordinates of a point of the 
system relatively to the rotating frame, 

\ OXy, OX^y 
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The equation of energy in general is at once seen to be 



dti. 



dt dt 












SO that if H be explicitly independent of the time, the forces Q,. be absent^ and P be 
a homogeneous quadratic function of x^, ,.., i„, 



"dt 



2f; 



§ 17. The simplicity of the fornmlation depends on the fact that the invariant 
factors of ^"'^a—xjr are linear. We have obtained this by assuming that the form ap^' 
only vanishes when every element of p is zero. But the invariant factors may be 
linear when this is not so, and the roots of the deterrninantal equation are not pure 
imaginaries. For instance, take Hill's equations for the motion of the moon, under 
certain limitations. 



d^x 



dv 



dt^ dt \r' 



0. 



y 



dt 



. ^ dx , K/M 

-T^i^^ "Tj -I 

(M 



0. 



/V 



r" 



riting 



li 



^n^oc^ 4- i (X -f nyf H- i ( Y — nx 



these are the same as 



X 

dt 



_„ 



oX ' dt 



3F 
dx ' 



y . d'F 

dt 



Q 



O 



Y' dt 



dv 



The so-called moon of no quadratures is obtained by variation from the solution 

expressed by 



lAj — {J » 



?/- 



= flcr. 



where o- is • given ■ by fx = Sn^o-^ ; this is a position of relative equilibrium. The 
matrix S- of the notation used above is zero ; the m.atrix a is 



^n^ 







n 



1 



n 







n 



An 







n 



i / 



In this case the quadratic form ap'^ is 



9n^p/ + { Pa + S^Ps)^ + h^Pi ^Pi)' 
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and vanishes when p^ =: p^ = 0, p2 " --2np^, But the roots of the determinantal 
equation 



are all different, and therefore the invariant factors are linear. 

The roots are 

^ =: + r4(28)'+lH V^ - ± m{(28)^-l}% 



of which only two are pure imaginaries ; thus not every disturbed orbit is periodic. 

I 18. We pass on now to give the details of the application of the general method 
above explained to the computation of some particular cases. 

A very simple case may be first given, merely as an example of the notation and 
method, since the results, once obtained, are easily verified. 

Take the equations 

= — ^ COS i^ + ^(l-f sin ^), 



Zi 



Zi 



dt 



These may be written 

d {x, y) __ 
d/t " 



dt 



^ (l — sin t)-\-y COS t. 



1 



\ 



1 



2/ 0, l\+2 



1, 



- COS t, sm t 
sin t^ cos t 




(«, y). 



or, say, 



d {x, y) 



a 



{xh-\-v){x, y), 



where 



u 



1. 

2 



0, l\, V 



2 



COS t^ sm A. 



have at once 



and therefore 



{2u)- 



~l, 



0, 



V-1, 




w 



W' .o 



i}(M)= l+t,,H-^«H-^7f/+... 



2! 



3! 



1 



sin t, cos t^ 



1 



1, 



r 



C08|-^+ / 0, 1\ W.\\t 



-1 



3! 



W+- 



--1, 



1 



cos -f^t, sm -i^t \ . 



sm 2''^3 



cos \t_ 

2 A 2 
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This gives 



Q-^ (t/,) =: /cos ^t, — sin ^t 



sin ijt^ cos ^t 



Wherefore 



Q ^ (u) . V ,Q (u) ™ |- /cos ^t^ — sin |-^ 



2 



sm ^?^, cos ^^ 




cos t, sin A 1} (u) 



sm ^. 



cos -1-^5 sin ^t\ Q (u) 



cos t 



sin ^^, cos -J^y 



i-/-i, 



0. 



1 



Denoting this by -Jo-, we find o-^ = 1, and hence 



— ch^t + arshi^t 



2^? 



Thus the solution is 



■u 



= /e-~^\ 0\. 



0. 




(^x, y) = Q{u-{-v) {x^, y^) = / cos ^t, sin -|A /e"'^^, 0\ {x^, y^), 



namely, 



sin ^t, cos ^t 



X — x^e"^'^ cos -^^-i-y^e^^ sin |-?^, 




2/ 



_ ^0^--|£ 



.o^it 



0, 



x'e"'^ sm ^t + i/e^ cos f^. 



The period of the coefficients in the original equation is Stt. The functions cos ^t, 
sin ^t have only the period 4x. To bring the result into the form given by the 
general theory we may write 

2* 



y = 



= - x»e-^'^+^" . ~. ((?'*- l) + ?/V<'+^'* . J (1 +e-**). 



2i 



the so-called characteristic exponents being 



±i (!+*)• 
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§ 19. We now consider cases of the equations 

these are derivable from the equation 

ft ^ 

df 
by taking 

leading to 

n 



As we wish particularly to illustrate the method of obtaining the characteristic 
exponents from the present point of view, we take first a case in which explicit terms 
in T arise early in the method of successive approximation. We take namely n ^ 1, 
and suppose 

^ — Xh + 2\ki cos 2t + 2\% cos 4^ + ... , 
4 

= Xh + XJciWi + X%2'^2 + - . ? 

where X is small, and w^ is used to denote i'' + f "''. 
Denoting ^yjr by ^, our differential equations are 

_ — ™ Cv \ .J-\.m JL /, 

(A/T 

where 

u = Y -—(p^ 0f\» 

The coefficients in these equations have period 2iri ; by what we have previously 
shown (§§ 14, 15), the solution is of the form 

(X, Y) = Vh/e-^\ \h-^ (K\ Y% 

0, e^' 



where P is a matrix whose elements have the period 27ri, h is a matrix of constants, 
and q is the constant which we particularly desire to find. As 

this corresponds to characteristic factors e^'^^+^^^^* for the original equation in t, whose 
coefficients have period tt. The quantity q is to be found by determming the terms in 
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T in the solution of the (X, Y) equations, and forming from this, after putting r — 0, 
a determinantal equation (§14), 

We are to calculate in turn Qii^ QiiQ'a^ &o,^ and arrange tlie result according to 
powers of X. First we have 



where 



a, 



C|, Clj/ 



fj> CvT^^ 



fr 



J " 



Ci = 



'<f> dr^ d'l 



m CtT I 
Jo 



thus, as (p is unaltered by changing the sign of r, bj can be obtained from c^ by- 
changing the sign of r, and similarly di from ai. This we denote by writing 



. . ^ aJL\_/4-JL 



and hence 



UUfil 



h 



<pl 



c 



If ^n 



of/ "I J. 



rv. ^> 




a 



15 '-' 1 



^h.5 ^ 







'?yC'Vv'.'?yC' / C|/0, C <, 



2? '^ 2 \ ? 



where 



Ci/( 






•JO 



p 



V- 



c\~\~ia\)dT, 



C2 = i'~^<l>{- % + ^^'i) dr, f/2 



</} {a 



^(-V,) dr, 



so that a^2 is obtained from ag by changing the sign of t- throughout, and similarly c'^ 
from Cg. In general, in passing from a term of Q (ii) involving r integrations to one 
involving (r+l) integrations, we shall have a law expressible by 






r+l 



<j>{-A, + ^G,)dT, 0,^1 = i^"V/.(-.A,+^0,)c?T, 



and the new term, like that from which it is derived, will be of the form 

-^r+1? ^ r + l 
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where A^^+i is derived from A^^^ by change of the sign of r, and similarly 0^^^^ from 

Thus, when in Q (u) we pick out the coefficient of r, as it occurs explicitly, 
independently of its occurrence in f, and in this coefficient put r — 0, we shall obtain 
a series of the form 



.71: 



Oti 



72. 



«%. 



where the first of these comes from Qu and involves terms in X and higher powers, 
the second comes from QuQu and involves terms in X^ and higher powers, and so on. 
And the equation for q will be of tlie form 



namely. 



061 + ^2+ ...—g, 

ri + y2+---> 



71 — 72—'* 



oti-as' 



-q 



(f = (ai + a2+...)^-(yi + y2+...)l 



Further, if the part of Qt^ ^hich is independent of explicit powers of r, consisting 
of elements which are polynomials in f, l"^ and periodic with period 27r^, be denoted 
by Pi, and similarly the periodic part of Qt^Qa be denoted by P^, &c., then the 
periodic matrix P above spoken of will be 

Proceeding to the computation, retain first only to terms in X. Then 



ai 



(^d 



T 



.XhT-^Xk.U-i'-'): 



C, 



Hence 



and q is given by 



Jo 



a^ 



r=\]l{t'-l)-VXh, 



™ --X/^, ji = — x^, 



't 



2 



In the case when the differential equation is that considered by Hill, this gives at 
once a very near approximation, as he remarks, being equivalent to his formula 



(Hill's ' Collect. Works,' I, p. 260). 
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If next we retain as far as X^ we have from 



= xh+\h{t'+o+^%{t'+n+^%{t'+n 



— a. 



d 



T 



xhT+\h{-r'+o+i'>^%{-t'+n+i'^%{-t'+n. 



Jo 



xh{t'-i)+\k(^-~ 



2 



■~'t)+X% 









i]+x%{^~+i 



r 



Hence 



r 



X/, (,+ 1 _f ) + x^, (- ^' +K-rf) + A^/% (- ^' +K- f ) + ^'^3 (- ^' + 



6/ 



rri 



Thus 



(p{-a, + ^c^) 






T — T^^ + ir)} 



This gives 

Jo 
^X^ih^r' + T+l-O + hk 



+hh{-U-'-rt'+ir'+^-^+it-'rn}. 



Similarly, 



-rt'-U-'-r + i + U-in 



Forming now —asf fcy we obtain 

A^ {A^ (-iT^-2T-3-T-(:+30 + M. iirt' + t' + r + h + M-S-U') 






/ 






DIFFEEENTIAL EQUATIONS OF ASTRONOMICAL INTEREST. 



m 



To find 



Jo Jo 



up to X^ it Is sufficient to take ^ = X [/^+^i (f^^+f)] ; so we obtain 

= X^ 1^3 ^_i^ «,2x— 3-— t|^+30 

S a'^S ^S +¥ 2^S 2*5 4^s ^S T-i-s / 



+ A^ii(— Jt^^™^ 



if) 









^3 = ^(-a2+fe)^T 

Jo 



+X%k{ 



( 



1 



^) ^-^+(^^^1) ^-1 + J^^3^|^_^|^^il^^^^_J^y^ 



+ 






4 



% 



6 






uiularly, 



<?3 = S^~"V (•-%+^<?2) <^T 







X%^[(|-t^+3t+6) i^""^ + (— Jt^ + St— 6)] 

+ X%^^j [(IV + T + I) r'-K-H(—|T«-iT^-4T-f) + (-T + |) ^] 

Picking ont- now the terms in t, putting therein f = 1^ and using the notation 
previonsly explained, we have, tip to X^^ 






'Xi 



If 



\^ ( - 2^^ + ^^i) + x^ iihh^ - hh). 



TOL. OOXVI.— A. 



2 B 



178 PROP. H. F. BAKER ON CERTAIH LINEAR 

and 

The product of these gives the value of q\ namely, 

This agrees with the value found above by a quite different method (§ 4). 
The matrix of coefficients of r, after ^ has been replaced by 1, is of the form 

A() =:/a, — y 

and its square is (a^— y^) times the matrix unity. The matrix Qq"^ (u) of § 14 is thus 

1 1 

0^^ (u) = 1 + Aq'm; + ig W + -~~: g^A^w^ + — g%* + . . . 

3 ! 4 ! 

or 

C + CtS, — yS X, 

s^ yS, G — aS^ 

where C = eh{qw), B = -sh (qw). From this it is easily seen that for the calculation of 
q the method we have followed is less laborious than to use the equation 

The differential equation from which we have started is, to terms in ^^^ if we 
suppose X ~ 1, 

-r-- 4- (l + 4h + 8ki cos 2^4-8Ecos it)x = 0. 

If we compare this with the form considered by Hill (' Coll. Works/ I.j pp. 246, 
268), we have, with his numerical values, 

h^ 0'039n 09848 99146, 
^^=^0-01426 10046 86726, 
k= 0*00009 58094 99389. 
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§20. Consider now briefly the case of the equations 

dr 4ni 



t-\ 1 



in which n = 2. We suppose 



An 



Xh + XkiWi + X^k2W2 + . . . , = ^5 say 5 



ttv-r^+r, 



As in the case ofn== 1, 



t^=/ -~^, f^^\. 



Qw 



r"<^, </> 



^ij ^ 1 \ J 



where 



a, = 



A 5 ^ 1. 



^(^T= -A/iT+x^,(r^-0+ix^'^3(r^-a, 



and, for n = 2, 



Ci = 







t'<pdT = xA(K-^-i)+x^i(K-'+r'-t)+x'^2(K-*-i-T). 



These give 



9^(-ai + rci) 



X%^ (r+i-l-r) +X%^i [(r-i) t'+ (t+2) f-K^-K'] 



As before 



QuQu 



+xV(-E-^+l-K+2r-|f). 



%J ^ 2 \ » 



where 



C2) €b 2, 



a^. 



Jss ■ 



+\%' {^^-^+%r+^-u+ r-if ) 



2 B 2 
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a.Tid 



c, 



ft 

Jo 






2/ 



+ X%l;i(— |-T^""'^ 



1 /^' 



■3 



rr^-sr^^-iT-K) 



+x^^i«(K- 



.4 S 



E-^+K-'+2T+i-| 



t). 



Picking out the coefficients of r in these, and putting therein ^ = 1, we have 



Ct J — — "**" A/tfj 



yi 



= -\'h 



2> 



a, = IX^^+iX^'^/, y, = -X%2-|X%Ai + 2X^/c/, 



and hence, to X'"^, 



T~ 



g2 :::^ ^^^ + ^^^2 _ ^^^ ^ ^^2 



\h 



|(F+|4^) 



9 



x^( m 



= x%^»^x%{F+p/). 



This agrees with a result previously found (| 6)^ but fails to" give the first term in 
q^ if Ji = 0. When this is so it, is necessary to take- account of the terms in \^. By- 
taking terms in X'^^'in %, Ci,-^-we only obtain terms. in ^(— c^i + ^^Ci) which involve X*. 
But the terms in X^ in ai^ Ci .which are written down give terms in X^ in ^ (— ^i + ^^c^i), 
which are 

and hence the additional terms in do 



x%^, (-K-'-f r^-ir+i+tr-K*) 






and the additional terms in c. 



/ 



i.™.i^~4 



9-2^ 



X'hk, ( — K 






3 . a^ I 1 



o 



L 
4 



+x^iA (Mr^^-Af-^^+f r^^+ii^^+rc 



,*— 1 



+-W 



rj-. 



(- 



si ' 



6 
T 



4\ 

§7* 



• In finding the terms in X^ in %, Og, it is sufficient to retain the terms X^ in % and Cg. 

This gives for ^ (—ag+C^Og), 



X¥?/ I 



2^ — T— 4i-4S • 



1"^^^} 



+x%^^arM-iT^~4-+A)"f+f(-iT«-ir-4)+f(-ir+^\a)+r(-iT+xy} 

;-x^/7/./{r^{|T+ii)~tri-fT-j.f+^(-~iT+Af)-^^r+r(-iT+i^)-ir} 

+x^/./{-l^«-^^(fT+f|)+|+r(fr-f|)^-wHf(2T-i)-|r}. 



DIFFERENTIAL EQUATIONS OF ASTRONOMICAL INTEREST. 



181 



For a. 



(p {—a2 + ^^C2)dT, this leads to 








12 
2^ 



2-T' 



3\ 
2/ 



+x%k,^-tHh+U)+it' 






-g-T 



+ S ( — 3T + -"9~) + S' 



1^'8 "^ S \ "3 T -t- -9--; 



T^ S \ 9^ r 2 7 / 



T- . 7 



6 



+ 3 6 



I2S 






The terms in Co — 



3 



f 2^ ^---^^^ 4.^2^2) (^^^ are similarly 



^'hHtniT^ + iT)+ir-ir'} 



+ \^h% 



^ '"^l 6T^+ — + T-S-8-)+3-f ^+^ ^. (,2 T^ + f T + -V") 






+ X%/^/ <^ f 



-4 



T 

6 



36 



)+ii '+s ("e^'7' + -;6'^-) + s (3'^ — '^^) 



"36 g-Ti-C> \^ — grT-t--3-; 



6^ 



+x«y^,M;Ay^■^+^^(l'T+li)-|^H^M-fT+f|) 

It is easy to see that the terms in X^ in a^i, Ci are respectively 

iA'U^3(r^-f) and X=;(;3{K-'-f+l), 
neither of which contains r. Thus up to X'^ we have, in the preceding notation 



C£i 



\K 



a,^ixHh'Hk^)~\'{hk, + ikJc,)., 



71 

72 



-X%, 



^3 = _ A%' _ %\%% - ^^--x^hk,' + ^i^k^', 73 = X'M + -'i-X%% - \^-X'hk' + '%^-X%\ 



ihus 



ai + as + «3 = -Xh + iX' {¥ + -P/) + X^{-V- %h% - 4f /i/t/ + \^'h'' - hk, - Uih), 
71 + 72 + 73 = X'{-h'--%hh + 2h'-h) + X' {h' + m%-^hk,' + -%'-h' + ihk,). 



182 PROF. H. F. BAKER ON CERTAIN LINEAR 

This gives 

(f = («! + 0^2 + ^3)^ -(71 + 72 + 73)^ 

- X%^-"X% (A^+fV) + X' (p' + ^^A'y^/-^#i'--y&/+4yfc//c,) 
as far as terms in X^ This result is for the equation 

wherein n ~ 2 and 



-\-{v?-\-\fr)x ^ 0, 



|-x/.+xy^,(r^+0+A^^2(r^+r)+... 

and agrees with the result previously found (§6) when in this last we replace 
h, hi, ^2 respectively by 2h, ik^, ik^, as is necessary, taking account of the difference 



when 



of notation for ^ in the two cases. By an independent investigation for the case 



we have found (above, p. 142), 

+ A' { {W ^W- h-^' + h' + ^i-h%' - 2Ai%2 + 2hji, - {k, - 2hJ } + ..., 

which, replacing h by hi + h^jX + h^X^^ arises from the preceding result. 
§ 21. Now consider the equations 

-^ (X, Y)^u (X, Y), 

CiT 

where 

and 9^ is not 1 or 2, but is an integer if u is sl periodic matrix. 
With 

^ = x/i+xA:i(r^+0+^'^2(r'+r) +... 

we have, retaining only to X^, 



a, — — 



<pdr= -Xhr + Xh ir'-OH-^'hit'-n 



c, 



n \7i + l n — 1 n—lj \n+2 n—2 n —4:/ 
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which lead to 



^(-ai + f^Cx) = X%^(t+ 



l~r' 



n 



r 



^X%hAt 



+ ^^K' 



T 



n{n-{-l) 



-4-f '7-4- -^ 



t?? + tl— 1\ ^ 



%-l 






n+1 



n{n'~-l) 

2ti 






n^—l n—1 W--1 



n n^—1 



o '2 -t ^ 



n 



so that 



%= ^{~ai+^Ci)clT 







X A ( -^T + — 



T r-1 



ti n 



n{n-\-l)/ \ n{n—l), 






^n-l 



2 



Pn 



t 



n+1 



n^—l n{n—l) n^—l n{n+l) 



J l/t a 



+ \% 



n 



t' + 



n 






2n ^ n^{n^—b) 

2% 



t 



9'%—%^ 



2n 



imilarly, 



2{n-\) {n-l){n'-l) {n+l){n'-l) 



w+1 



C2= rv(-ai+rci)<^T 



\ 



21,2 



\% 









,2 






2Z. 2 



+ X')fci 



n 



(^+l) (n + 2) 



&-.n-~2__ -^ f -n 



n 



n'-l {n-l){n-2) 



r 



•n + 2 



+ ' - :■ S ^ ; — 7 f + 



B^— 1 fi-^— 1 ti^— 4 



Thus we have, so far as terms in X^, 



eti 



Ail. 



ri = o, 



a. 






71- 



n^— 1' 



72 



2X%^ AX^hjn 
n n^— 1 ' 
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and so, to this approximation, 



n n-~l 

The characteristic factor is then e^^^'+^^^^^ the differential equation being 



dt' 



+ {n^ + inXh + SnXJci cos 2t -f SnX^ks cos U-i- ...)x ~ 0. 



Thus q is always real, when X is small enough, provided ki is not zero, even if h be 
zero. The result agrees with that found in § 6 for n ^- 3^ if allowance be made for 
the change of notation. 

[^Decemhef' i, 1915. — -Consider the differential equation differing from that just 

preceding only by the substitution of H for \h in the term inXh of the coefficient of 

X, where H is supposed to be of the form X A^ + X%2 + ^^^3 + — Th© computation of q^ 

proceeds then exactly as before. The formulae for cci + 062 + ^3, 71 + 72 + 73, given above, 

p. 178, substituting H for \h^ show that, for n = 1, q^ is then of the form 

(H — cti) (H — ^2) Q, wherein Q is a power series in H, X^i, X^^2? •••? reducing to 1 when 

H = 0, X ™ 0, and 

a^=: - k^X - l^i^X^ + {ik{^ - ^1^2) X'V . . . , 

^2 — kiX — "l^/X^ — {^kj^ --kJc2)X^ -{-.... 

The value of q^ is positive, and the motion represented by the differential equation 

is stable, so long as H does not lie between these values. Similarly for n ™ 2, from 

the formulae at the bottom of p. 181, the range in which q^ is negative is when H lies 

between 

-(Pi'"^2)X' and (YA'-^2)a^ 

these being accurate as far as X'l Unless ^k^ < ^2 < ^i~ki\ these limits are of 
opposite sign, and include H — 0. This is the result given on p. 142 (save for a slight 
difference of notation). For ^ = 3, an analogous computation shows that q^ is positive 
except when H is between 

3^/X^^PX^ and f^,V + PX^ 

where 

V ^W---3kA-^h, 

and this range does not include H = unless ki ™ 0. It would appear, from the 
formula above (p. 184), that the corresponding interval for greater integer values of ^ 
is between two quantities of the forms 

^"^ k;'x'^-vx\ -l^^iV+Qx", 



n'-l ' ' n'-l 
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Mr. E. Lindsay Inoe, of Trinity College, Cambridge, following up the method of 
his paper referred to above (footnote, p. 134), has calculated numerical results for the 
case when hi^ ks, ... have the values considered by G. W. Hill._ 

PAT^T TTT 

§ 22. I desire to add to the foregoing some very incomplete remarks in regard to a 
generalisation of which the work appears to be capable. The most important general 
result obtained is that when ti is a periodic matrix, the matrix Q (it) can be expressed 
as a periodic matrix P multiplied into a matrix involving quantities of the form e^"". 
One direction in which this result can be amplified is by extending the assumption 
we have made that the matrix Qq^{u) has liDoar invariant factors. It is well enough 
understood what is the character of the modifications thereby introduced. A more 
important generalisation appears to be that the factorisation of the matrix Q{u) 
does not in fact require that u be a periodic matrix. As an indication of the theorem 
consider an equation 

in which the constants k, \ jn are such that /c+X + m = 0, but the ratio of two of them 
at least is irrational. For example, we might have /c='\/2 + l, ,u=— v^2 + l, 
/x = — 2. Then, assuming that there exists no identity of the form 

a/c + /3X + y/Ai2cr = 0, 

in which a, /3, y are positive integers, the equation would seem to have a solution of 
the form 

where X is a series of positive and negative integral powers of e*"*, e'^*, e^^\ e''^\ which 
may be arranged as a power series in a, h, c, and q is a series of the form 

q = (r+A]a5c + A2a^6V+..., 

in which Aj, Ag, ... are constants. The difierential equation has not periodic 
coefficients. 

In a paper already far too often referred to, ^ Proc. Lond. Math. Soc.,' XXXV., 
1902, p. 353 et seq.^ replacing the variable there called t by e" or ^, it is shown 
(p. 365) for the equation system 

— = (A + fV)^, == ux, say, 

in which A is a matrix of constants, and Y a series of positive integral powers of f, 
that there is a factorisation of the matrix 0(fi), in the form P0(9i)y, where P is a 
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matrix whose elements are power series in f, and {<p) is calculated in regard to 
from a matrix 



a— '^s 



while y is a matrix of constants. Here 6i, 62^ . . . , c^o, c^^^ . . . depend solely on the 
invariant factors of the matrix A. 

This result is obtained from the form of the matrix u as expressible by powers of f, 
without reference to the question of periodicity. It would seem that the argument 
there employed is capable of modificatioUj the integrations being performed in regard 
to T (which is log ^ of the paper referred to)^ so as to lead to the general theorem 
here contemplated. 
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